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Topological insulators and topological superconductors (TSC) display various topological phases 
that are characterized by different Chern numbers or by gapless edge states. In this work we 
show that various quantum information methods such as the von Neumann entropy, entanglement 
spectrum, fidelity and hdelity spectrum, may be used to detect and distinguish topological phases 
and their transitions. As an example we consider a two-dimensional p-wave superconductor, with 
Rashba spin-orbit coupling and a Zeeman term. The nature of the phases and their changes are 
clarified by the eigenvectors of the /c-space reduced density matrix. We show that in the topologically 
non-trivial phases the highest weight eigenvector is fully aligned with the triplet pairing state. A 
signature of the various phase transitions between two points on the parameter space is encoded in 
the /c-space fidelity operator. 

PACS numbers: 03.65.Ud, 03.67.Mn, 74.40.Kb, 03.65.Vf 


I. INTRODUCTION 

Topological phases are non-trivial phases of matter 
characterized by global entanglement and correlations. 
Various examples have received increased attention such 
as the traditional quantum Hall phases, spin liquids, 
topological insulators and topological superconductors. 
Due to their nature, topological phases are robust to lo¬ 
cal perturbations and have received wide attention as 
possible elements for error-free quantum computation. 

In two dimensions, time reversal breaking Z topolog¬ 
ical insulators exhibit a charge Hall conductivity that 
is quantized and proportional to the Chern number of 
the occupied band^. Time reversal invariant Z 2 topo¬ 
logical insulators have also been proposed and found 
experimentall}BM^. Such nontrivial topological phases 
are also characterized by the presence of gapless edge 
mode^^. 

Superconductivity with non-trivial topology may also 
be obtained^. It can be due to the pairing symmetry, 
as is the case of p-wave SCJ^. In semiconductors with 
Rashba SO coupling it arises when s-wave superconduc¬ 
tivity is induced and a Zeeman term is adde(P®2. In case 
the normal phase is already topologically non-trivial, a 
TSC may be obtained if g-wav e superconductivity is in¬ 
duced by proximity effecfP^^. Topological phases and 
their phase transitions in topological superconductors 
may be detected through the existence of zero bias peaks 
in tunneling spectroscopy experiment's®^ due to Ma- 
jorana end states, including spatially resolved peakJSl^ 
or through anomalous Fraunhofer patterns or fractional 
Josephson effectJSMU. Multiple Andreev reflection cur¬ 
rent in voltage-biased junctions has also been proposed 
as a signature of topological ordei!^, in particular, in the 
derivative of the current and in the zero-bias conductance 
of nodal noncentrosymmetric SC^Sl^ xhe imaginary part 
of axy at finite frequency may also be used to signal topo¬ 
logical phased, as well as the Hall conductivity and its 
derivativeP^. 


On the other hand, the interplay between quantum in¬ 
formation and condensed matter physics has been exten¬ 
sively considered, such as the use of entanglemenfP^in the 
study of zero-temperature quantum phase transitions^. 
This interrelation has been explored in the reanalysis 
of several non-trivial exactly solvable models, using dif¬ 
ferent information measures to better understand the 
underlying physicS^®^. For a bipartite system, be¬ 
sides the von Neumann entropy and related quantitieS^, 
other information measures like the concurrencS^, 
the mutual informat iorP^®^, the negativit}02, or the 
Meyer- Wallacl P^ an d the generalized global entangle¬ 
ment measureS^^®! have been considered. The effect 
of the quantum statistics has been analyzed for free 
electrons^ and bosonS^, and for electrons in a BCS 
superconductoJ^. 

The distinguishability between states through the fi¬ 
delity, has also been used as a possible criterion to de¬ 
tect quantum phase transitions. By its own nature, fi¬ 
delity between pure ground states signals a change of 
state as one approaches a quantum phase transitioiP^®^. 
The fidelity between mixed states has als o bee n used 
as a signature of quantum phase transition^^^^^ and to 
distinguish between different states of matter at finite 
temperature^®®^. 

As argued in Ref. [64l a more detailed information 
about a mixed state may be obtained if the entangle¬ 
ment spectrum is analyzed. Considering reduced density 
matrices, where part of the degrees of freedom are inte¬ 
grated over, it was shown, in the context of the quan¬ 
tum Hall effecfP^^ and in the context of coupled spin 
chainJ®!, that the ground partitioning state entanglement 
spectrum of a subsystem A contains information about 
excited energy states, of the frontier of the subsystem 
A with the complementary subsystem B. Other parti¬ 
tioning of the system have been proposed that lead to 
further informatioi]®l. Considering a partitioning in mo¬ 
mentum space, it was shown that information about en¬ 
ergy excitations of a single Heisenberg chain is contained 
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in the groundstate wave function, through the entangle¬ 
ment spectrunP^. 

One may also considethe spectrum of the fidelity 
operator J-{pi^p 2 ) between two density matrices, pi and 
p 2 . Its set of eigenvalues which is denoted fidelity op¬ 
erator speetrum^ and — In/^, called the fidelity speetrum^ 
provide more information as compared to the fidelity (its 
trace). This parallels the extra information provided by 
the entanglement spectrunP^, as compared to the von 
Neumann entropy. In the case of two equal mixed states, 
the operator T has a set of eigenvalues, = A^, such 
that — In reduces to the entanglement spectrum. 

Quantum information methods have also been used 
to detect the more elusive topological phases and tran¬ 
sitions. Kosterlitz-Thouless transitions were success¬ 
fully detected cal culat ing the fidelity susceptibility of the 
XXZ spin chaiii^^^. A study of the one-dimensional 
asymmetric Hubbard model showed that, different scal¬ 
ing regimes of the fidelity susceptibility may be used to 
distinguish the two phased. Other systems that display 
topological phases were also studied, like spin-1/2 par¬ 
ticles on a toru^^, the toric code model and the quan¬ 
tum eight-vertex modeP^ and the spin honeycomb Ki- 
taev modeP^. In the case of the Hubbard model it was 
shown that the fidelity metric satisfies an hyper scalin^^ 
and in the Kitaev honeycomb model it was found a differ¬ 
ent temperature scaling behavior in the different phase^^ 
and a divergent fidelity per lattice site, showing that a lo¬ 
cal measure is also able to detect the global entanglement 
regim^^, as in other system^^^H^. 

In this work we will be interested in applying various 
quantum information measures to study the topological 
phases and phase transitions in a topological supercon¬ 
ductor. In particular, we will consider the entropy, en¬ 
tanglement spectrum, fidelity and fidelity spectrum, fo¬ 
cusing on global properties. We consider a momentum 
space partitioning, since the Hamiltonian factorizes. In 
general, the groundstate is a singlet and the entropy van¬ 
ishes. At the transitions between the various phases the 
spectrum is gapless and degenerate leading to a finite en¬ 
tropy. The opposite result is obtained for the fidelity be¬ 
tween two density matrices, defined close by in parameter 
space; when both density matrices are defined in the same 
phase the fidelity approaches one and when the transi¬ 
tion is reached the fidelity has a minimum. However, 
the momentum space contributions to both the entropy 
and fidelity, or the entanglement spectrum and fidelity 
spectrum contain a richer structure. Also, the eigenvec¬ 
tors of the reduced density matrix reveal the nature of 
the state in a given phase and, in particular, the change 
across a topological transition from a trivial to a non¬ 
trivial phase. A recent use of the fidelity and fidelity 
susceptibility was carried out to study a one-dimensional 
spinless superconductor with end Major ana fermion 
Here we consider a two-dimensional superconductor with 
a parity breaking Rashba spin-orbit term that, in general, 
mixes singlet and triplet pairingJ^. Also, a Zeemann 
term is added leading to a rich phase diagranP^ with 



FIG. 1: (Color online) Topological phases and their Chern 
(C) numbers as a function of chemical potential and mag¬ 
netization. The transitions occur at three sets of momenta 
k = (0, 0), k = (tt, 0), k = (tt, tt). 


different topological phases, that can be characterized 
by the Cher n number, gapless edge states and winding 
numbeil^^^. 

In section H we present the model studied and the en¬ 
tities used to detect the quantum transitions. In section 


HI we present the calculation method and our results 


for the entropy. In section |TV| we present results for the 
entropy spectrum and the analysis of the band struc¬ 
ture of the system. Also, we show how the eingenvectors 
of the system can be useful for a better understanding 
of the transitions. The eigenvalues and eigenvectors of 
k-subspace reduced density matrices are considered in 
subsection |IV A[ In section [V| we present the calculation 
method and our results for the fidelity, and show how it 
clearly detects the quantum transitions. In section [Vl| we 
present our main conclusions. Finally, in the appendices 
we show in detail the results for the eigenvector of the 
largest eigenvalue of the density matrix and the reduced 
density matrix. 


II. TOPOLOGICAL SUPERCONDUCTOR 


We consider a two-dimensional triplet superconductor 
with p-wave symmetry. This model was studied in Refs. 
I25l^ We write the Hamiltonian as 




where (4,c_fe) = (4t, C-fet, C-fe^) and 




( 1 ) 


Ho = efcCTo - + Hr ■ (2) 

Here, = —2t{coskx + cosky) — Sp is the kinetic part, 
t denotes the hopping parameter set in the following as 















3 


the energy scale, ei? is the chemical potential, /c is a wave 
vector in the xy plane, and we have taken the lattice con¬ 
stant to be unity. Furthermore, Mz is the Zeeman split¬ 
ting term responsible for the magnetization, in energy 
units. The Rashba spin-orbit term is written as 

Hr = s • cr = a (sin kyCTx — sin kxCJy) , (3) 


where a is measured in the energy units and s = 
a{sm ky, — sin kx^O). The matrices (Tx^dy^az are the 
Pauli matrices acting on the spin sector, and do is the 
2x2 identity. 

The pairing matrix reads 


A = i{d- a As) dy 


dx T idy dz T Ag \ /^\ 

dz Ag dx T idy J 


We consider a situation where the spin-orbit coupling is 
such that the pairing is aligned^ along the spin-orbit 
vector s as d = {dx^dy^dz) = {d/a)s and d is a scale 
parameter. This is a situation expected if th e spin-orbit 
is strong (other cases were considered 

The energy eigenvalues and eigenfunction may be ob¬ 
tained solving the Bogoliubov-de Gennes equations 


f Ho{k) A{k) (u, 

V At(fe) -Hu-k) J G- y '’’"V 


( 5 ) 


The 4-component spinor can be written as 

/ Un{k^^) ^ 

f A _ 

\^n J ~ Vni-k,^) 

\ Vn{-k,i) / 


The energy eigenvalues can be written as 


^k,n 


ai\Jzi+ 2q!2v^, 


where <^1^2 = =tl and 

Z\ = d*d + s- s + 6^ + -|- Ag 

and 

^2 = ^ ^ ^)zf 

■T Ag (^Agi^di • d + + 2efcd • s) . 


(6) 

( 7 ) 

( 8 ) 

(9) 


The gap of the lowest band closes when zi = 2yT2. 

Since a topological transition may occur when a gap 
closes, the gapless points may indicate the presence of 
topological transitions. One way to characterize various 
topological phases is through the Chern number, obtain¬ 
able as a n int egral over the Brillouin zone of the Berry 
cnrvatur cj^^ l ^^l The locations of the gapless points may 
also be detected by the Berry curvature of the bands. 
When two bands become close there is a peak in the 
Berry curvature. It can be shown that for the lowest 
band, there are sharp peaks at the time-reversal momenta 
(as defined in Ref. 15^ /c = (0,0), Aj = (7r,0), (0,7r) and 


k = (7r,7r). The second band shows peaks at the vari¬ 
ous characteristic momenta, as discussed above for the 
entropy. Summing over th e occ upied bands the Chern 
number has been calculatecP^^^. The results in the pa¬ 
rameter space are shown in Fig. [fusing the typical pa¬ 
rameters t = 1, a = 0.6, d = 0.6 and = 0.1. From 
now on all figures use the same set of parameters, unless 
stated otherwise. 

Non-trivial topological order for non-interacting 
Hamiltonians can be related with the presence or absence 
of time-reversal symmet ry (T RS), particle-hole symme¬ 
try and chiral symmetr}®^. In Bogoliubov-de Gennes 
systems, particle-hole symmetry is always present, and 
TRS is determinant to the nature of possible topological 
phases in two dimensions. The superconductor we con¬ 
sider here is time-reversal invariant if the Zeeman term is 
absent. The system then belongs to the symmetry class 
Dill where the topological invariant is a Z 2 index (the 
TRS operator T is such that = — 1). If the Zeeman 
term is finite, TRS is broken and the system belongs to 
the symmetry class D (the TRS operator T is such that 
7^2 _ topological invariant that characterizes 

this phase is the first Chern number (7, and the system 
is said to be a Z topological superconductor. 

Due to the bulk-edge correspondence, complementary 
information on the topological phases and transitions 
may be obtained analyzing the edge states. If Mz = 0 
and the pairing is s-wave, the system is in a topologically 
trivial phase: there is only the bulk gap and no gapless 
(edge) states. In the case of p-wave or when there is a 
mixture of s- and p-wave components, and the ampli¬ 
tude of the p-wave pairing is larger than the correspond¬ 
ing amplitude of the s-wave case, there are two counter- 
propagating edge modes that give opposite contributions 
to the total Chern number, (7 = 0 (Z 2 phase]P^. 

As the magnetization is turned on, TRS is broken. For 
small magnetization, the superconductor is in a trivial 
phase with Chern number (7 = 0. A finite magnetization 
is then necessary to cause a topological phase transition 
to a phase with non-zero Chern number^. This happens 
both for the p-wave case and the s-wave case. The se¬ 
quence of Chern numbers is clearly correlated with the 
number of pairs of edge states as shown in Ref. [85l It is 
interesting to note that even though the system is in a 
(7 = 0 phase the number of edge states is two; the same 
as that in the parent Z 2 phase, when Mz = The 

presence of edge modes induced by bulk topology can 
also be shown using dimensional reduction and thereby 
calculating the winding numbeJ^. The calculation of 
the winding number gives the number of gapless edge 
modes both when the Chern number vanishes and when 
the Chern number is finit^^. 


III. VON NEUMANN ENTROPY 


Since the Hamiltonian is separable in momentum 
space, the density matrix operator for a momentum k 
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FIG. 2: (Color online) Total entropy as a function of chemical 
potential and magnetization, for the typical parameters t = 1, 
a = 0.6, d = 0.6, As = 0.1 and temperature T = 10“^. 


may be defined as usual as 


Pk = 


Q-PHk 


( 10 ) 


with H = (1/2) Ff/c, defined from equation (jlj. In 
the diagonal basis it is written as 


Pk = {n\pk\n) = 


-p(^n\Hk\n) 

Zk 


( 11 ) 


It is convenient to introduce a basis representation for 
the density matrix in terms of the occupation numbers 
for a given momentum (and its symmetric) and the two 
spin projections. The eigenvalues of the density matrix 
are obtained if we diagonalize the Hamiltonian in the 
same basis. We consider the representation 

Hk — /c^ I -^/c I'^k^'^—k^'^ki'^—ki^ (1^) 

The diagonalization of the Hamiltonian matrix in this 
enlarged basis is written as 


HkQ k,n — ^k,nQ k,r 


n = 1,..., 16 


(13) 


note that n here is just an index number and should not 
be confused with the occupation number of equation (12). 
In the same basis the density matrix may be written 


as 


Pk 


Q-/3Hk 


(14) 


Therefore the eigenvalues of the density matrix may be 
written as pkQk,n = ^k,nQk,n where 


A 




/c,n 


Ee 




(15) 


The entropy for each momentum can be obtained using 
that 


16 


^k — ^ ^ ■^k,n 


(16) 


n=l 


and the total entropy is obtained summing over the Bril- 
louin zone 


S = '£Sk. (17) 

k 



K 


FIG. 3: (Color online) Bands of Hk for Mz =0.1 and = 
—4, 0,4, respectively. 

The results for the entropy as a function of chemical 
potential and magnetization (Zeeman term) are shown 
in Fig. for a temperature parameter T = 10“^, where 
P = 1/T {ks = At very low temperature the en¬ 
tropy is of the type S' = Infl, where ft is the degeneracy 
of the lowest state. In general, the lowest state is non¬ 
degenerate and the entropy vanishes. At the transition 
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^k=(0,0) 


^k=(7t,0) 



FIG. 4: (Color online) Entropy as a function of chemical 
potential and magnetization for k = (0,0), k = (7r,0),k = 
(7r,7r), respectively. 



points the system becomes gapless and, at some points 
in momentum space, the lowest state becomes degener¬ 
ate and the entropy is finite. Integrating over momentum 
we sum over the various gapless points. As can be seen 
comparing Fig. [^and Fig. [^the entropy tracks very well 
the topological phase transitions. Except at the transi¬ 
tions, the entropy is featureless and does not distinguish 
between the various topological phases. 

Further insight is obtained if we analise the entropy 
Sk as a function of momentum, in order to understand 
the origin of the peaks in Fig. At Mz =0.1 there are 
three peaks in Fig. |^at chemical potentials ep = —4, 0,4. 
Those are special points, degenerate at absolute zero 
temperature, with values ln(4), 21n(4) and ln(4), respec¬ 
tively. The fact that those points occur at Mz = 0.1 has 
to do with the choice of = 0.1 (for = 0 those peaks 
occur at Mz = 0). These peaks are due to singularities at 
momenta k = (0,0), Aj = (7r,0),A5 = (7r,7r), respectively, 
and their equivalent points in the Brillouin zone. At very 
low temperature the peaks get very sharp. Therefore, a 
small finite temperature is used for better visualization. 

The energy spectrum degenerate gapless points are also 
seen in the bands of the matrix representation of the 
hamiltonian These are shown in Fig. where we 
have considered the same special points, Mz = 0.1, in 
the phase diagram. Their relevance also extends beyond 
those particular chemical potential values. As shown in 
Fig. I^the momenta responsible for the various transition 
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FIG. 5: (Color online) Eigenvector components for the high¬ 
est eigenvalue of the reduced density matrix pk for Mz =0.5 
as a function of the chemical potential ef for different momen¬ 
tum values k = (0,0), k = (tt, 0), k = (tt, tt). The components 
that are not shown are nearly zero. 
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FIG. 6: (Golor online) Eigenvector components for the high¬ 
est eigenvalue of the reduced density matrix p/e for = — 1 
as a function of the magnetization for different momentum 
values k = (0,0), k = (7r,0),k = ( tt , tt ). The components 
that are not shown are nearly zero. 


lines is the same set. The lines that emerge from sp = —4 
are associated with k = (0, 0), the lines that emerge from 
= 0 are associated with k = (7r,0) and those from 
= 4 with k = ( tt , tt ). 
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FIG. 7: (Color online) Eigenvector components of largest 
eigenvalue of reduced density matrix pk for arbitrary momen¬ 
tum values for Mz — 0.5 as a function of chemical potential 
(top panel) and for — 1 as a function of magnetization 

(lower panel). The components that are not shown are nearly 
zero. 


IV. ENTANGLEMENT SPECTRUM AND 
EIGENVECTORS 


The entanglement spectrum of a bipartite system with 
a real space partition, provides interesting information 
about the states along the interface between the two sub¬ 
systems, as discussed in the introduction. The momen¬ 
tum space partition carried out in this work, originating 
in the Fourier partial diagonalization of the Hamiltonian, 
provides information about the spectrum in each momen¬ 
tum separately. Even though this partition is trivial, 
a detailed analysis reveals interesting information about 
the entanglement of the system, through the entangle¬ 
ment spectrum and the eigenvectors of the reduced den¬ 
sity matrix. The entanglement spectrum can be seen 
from Fig. since it is simply the spectrum of the hamil- 
tonian matrix once both are diagonal in the same 
basis states. At low temperatures, this translates to a 
spectrum of the /c-space density matrix where the low¬ 
est band has an eigenvalue close to unity and all others 
close to zero (except at the degenerate points where the 
lowest state has an eigenvalue which is the inverse of the 
degeneracy). 


The eigenvectors provide additional information about 
the phases of the system, and in particular provide inter¬ 
esting information about the transitions. In general, each 
eigenvector is written as a linear combination of the 16 
basis states used, It turns out that in 

most cases, out of the 16 coefficients, only a few are non¬ 
zero, for each eigenvector. Their physical interpretation 
is in several cases clear, and it is possible to understand 
the physical content of the dominant eigenstates, and 
what changes as a transition occurs. 

The dominant states are: i) the states |1100) and 
10011) corresponding to the spin triplet pairing, ii) the 
state 10000) corresponding to an empty state in momen¬ 
tum and spin spaces, hi) the state |1111) corresponding 
to a basis state fully occupied with 4 electrons (for a 
given momentum and its symmetric and both spin com¬ 
ponents), iv) the states |1001), |0110) corresponding to 
the spin singlet pairing. Other states are also present 
but their coefficients are typically quite small. 


In Figs. mzi we show results for the absolute value 
of the coefficients of the most important contributions to 
the eigenvector, corresponding to the highest eigenvalue 
of the k- density matrix (corresponding to the lowest 
band of the hamiltonian matrix). In Fig. [^the magne¬ 
tization is fixed, the chemical potential is changing and 
each panel corresponds to the three singular momenta 
values. The set of chemical potential values are those 
corresponding to the horizontal dashed line in Fig. In 
Fig. 1^ the chemical potential is fixed at ep = —1 and 
the magnetization is varied tracing the vertical dashed 
line of Fig. In Fig. we consider similar plots but 
for arbitrary momenta values not corresponding to any 
of the singular momenta. 

Along the dashed horizontal line the Chern number 
changes from C = 0^C=l^C = 0-^C=—2^ 
C = 0^C = 1^C = 0. The symmetry in the se¬ 
quence of the Chern numbers around half-filling {sp = 0) 
is however not directly manifested, neither in the results 
for the coefficients, nor for the momenta where the tran¬ 
sitions are detected. But the symmetry comes out if 
we group the coefficients by similar spin pairing states, 
as shown in the appendix. The transitions crossed by 
the dashed line are detected considering k = (0,0),/c = 
(0,0),/c = (7r,0),/c = (7r,0),/c = (7r,7r),/c = (7r,7r), re¬ 
spectively. This is understood from the results of Fig. 

Bl 


The most important contribution of the coefficients 
analysis is that they reveal the nature of the states in each 
phase. An individual special k analysis shows that, for 
small chemical potential values, the empty state |0000) 
has the highest coefficient. As a transition occurs, the 
spin triplet basis state |1100) becomes the dominant until 
the next transition occurs, when the fully occupied state 
I fill) becomes the dominant one. This sequence hap¬ 
pens at different chemical potentials: those correspond¬ 
ing to the specific transition line for which, at given mo¬ 
mentum value, the lowest band becomes gapless. Close 
to the various transitions the spin singlet states |1001) 
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and 10110) also contribute. 

Similar results are obtained crossing the transitions by 
the vertical dashed line, shown in Fig. where the se¬ 
quence of transitions is now observed as a function of 
varying magnetization, at fixed chemical potential. How¬ 
ever, it is shown that, after a given transition, the spin 
triplet state becomes the dominant one and prevails un¬ 
til large values of the magnetization. Also, for low Mz 
and before the transition, depending on the momenta 
considered, the dominant state is either the empty state 
or the fully occupied state. It turns out that, counter¬ 
intuitively, for chemical potential sp = —^ the dominant 
state is the fully occupied state for momenta close to the 
origin of the Brillouin zone, while the empty state is the 
dominant contribution close to the edges of the Brillouin 
zone. 

The sharpness of the contribution of the various coef¬ 
ficients is lost, as shown in Fig. if the momenta are 
not the singular points. In this case no transitions are 
observed but their effect is still seen in broad features. 

As Figs. mm indicate, the topological phases with fi¬ 
nite Chern number, are associated with a large coefficient 
of the triplet contribution. In Fig. |^we show the mo¬ 
mentum dependence of the absolute value of the triplet 
pairing coefficient for Mz = 0 (left column) and Mz = 1 
(right column). Note that for Mz = 0 the system is in a 
Z 2 phase in the region where \sf\ <4. For finite magne¬ 
tization the system is in a Z phase with either finite or 
vanishing Chern number. For both values of the magne¬ 
tization, there is a clear correlation between the absolute 
value of the triplet coefficient and the location of the 
Fermi surface which expands as the chemical potential is 
increased. At small chemical potential the Fermi surface 
is centered around zero momentum. Increasing the chem¬ 
ical potential the Fermi surface will approach the point 
k = (tt, 0) and equivalent, and as the chemical poten¬ 
tial increases further the Fermi surface will approach the 
point k = ( tt , tt ) and equivalents. The amplitude of the 
coefficient is larger in the case of Mz = 1. For Mz = 0 the 
two spin orientations of the triplet state are degenerate 
and it can be seen that the amplitude is fairly constant 
along the Fermi surface. When the magnetization is fi¬ 
nite, there is a larger amplitude in the neighborhood of 
the singular points previously found for the entanglement 
entropy and fidelity. In this case, near these points, and 
depending on the chemical potential, the coefficient is 
close to saturation, which implies that the triplet state is 
dominant over all other basis states. 

In Appendix A further detailed results for the eigen¬ 
vector structure are presented. The results show that the 
topological phases may be identified by the saturation of 
the spin triplet contribution. The empty and doubly oc¬ 
cupied states usually dominate over the Brillouin zone, 
except near the Fermi surface, where the main contribu¬ 
tion comes from the triplet and singlet pairing states. 



FIG. 8: (Color online) Absolute value of the coefficient of 
the spin triplet state as a function of momentum. In the left 
column Mz = 0 and in the right column Mz — 1. The rows 
are for = —4, —2, 0, 2,4, respectively from top to bottom. 


A. /c-subspace reduced density matrix 

Since the Hamiltonian is additive in momentum space 
the density matrix is factorizable leading to a density 
matrix as expressed in eq. (11). We have diagonalized 
using the number occupation basis shown in the same 
equation. This leads to a 16 x 16 matrix both for the 
Hamiltonian and the density matrix (this density ma¬ 
trix may be understood as a reduced density matrix by 
integrating over all other momenta). Even though the 
many-body system is highly entangled in terms of wave 
functions (Slater determinant), at zero temperature this 
entanglement is hidden by the occupation number repre¬ 
sentation. We may however, consider a subspace within 
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FIG. 10: (Color online) Entropy as a function of momentum 
for Mz = 4, = 2 for the reduced density matrices p and 

p", respectively. 
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FIG. 9: (Color online) a) Lowest energy band, b) largest 
eigenvalue of reduced density matrix and c) entropy as a func¬ 
tion of momentum for Mz = 2,£i? = 0 (C = —2). 


each momentum value. 

The basis states lead to a deusity matrix expressed as 



M=2.8 - M=3.1 

M=2.9 - M=32 

M=M=3.0 



p(^/c,t5 ^ —/c,t5 ^ —/c,t5 ^ —/c,t’ ^ —/c,t) 

where each uumber operator takes the values 0,1. The 
diagoualizatiou of the Hamiltouiau was obtaiued through 
a uuitary trausformatiou iu the way 

Hd = QHQ-^ (18) 


This trausformatiou also diagoualizes the deusity matrix, 
which allows us to calculate the diagoual form of it, pd^ 
through eq. (15). Therefore, haviug diagoualized the 
Hamiltouiau we cau obtaiu the deusity matrix 

p = Q~^PdQ (19) 


FIG. 11: (Color online) Entropy as a function of kx for a) 
£f = and ky = 0 and b) for = 1 and ky = tv for values 
of Mz around the topological transitions. 


We may now construct a /c-subspace density matrix 
with dimension 4 x 4 in a reduced space by tracing out 
part of the degrees of freedom. If we define as part A 
the k momenta (with up and down spin) and a part 
B with momenta —k (with up and down spin), the re¬ 
duced density matrix is obtained tracing out the degrees 
of freedom, for instance, of part B, and we are left with 
a density matrix for part A. This choice gives the en- 
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tanglement between A and B which is the entanglement 
between electrons with momentum k and momentum —k, 
characteristic of superconductivity, including the singlet 
and the triplet pairings. 

Let us call p' the reduced density matrix of dimension 
4x4. It can be obtained element by element as 

p'(nfe,t, nk,i-, 

( 20 ) 

This density matrix gives information on the entangle¬ 
ment between electrons of opposite momenta. 

Alternatively, we may also define 

,4r 5^ — fc ,4' 

^ —/C,i"5 ^ —/c, 4.5 ^ —/C,t’ ^ —^d) 

( 21 ) 

and now, this density matrix gives information on the 
entanglement between electrons of opposite spins. 

The diagonalization of these reduced density matrices 
leads to 4 eigenvalues. Typically the two smaller eigen¬ 
values take very small values and, since they are nor¬ 
malized to one, the second largest eigenvalue is basically 
the difference to one of the largest eigenvalue. As shown 
befor^^ the largest eigenvalue takes values close to unity 
in most of the Brillouin zone, except near the Fermi sur¬ 
face, where it decreases in value. This is shown in Fig. 
in the first panel we plot the lowest positive energy band 
obtained from the solution of the Hamiltonian, eq. 0 
and in the second panel we plot the largest eigenvalue of 
the reduced density matrix p'. Further information may 
be obtained calculating the entropy of this /c-subspace 
reduced density matrix, defined as 

4 

^k = ~ ^ ^k,n l^(^/c,n) (22) 

n=l 

The results for the entropy are shown in the third panel 
of Fig. I^for the same set of parameters, = 2,5^ = 
0. The decrease of the largest eigenvalue is matched by 
the increase of the entropy, which therefore locates the 
region of the transition in momentum space between the 
occupied states and the unoccupied states defining the 
Fermi surface. 

In Fig. we compare the entropy between the re¬ 
duced density matrices p' and p" at another point in the 
phase diagram, = 4,5^ = 2 {C = —1). The entropy 
is clearly larger in the case of the entanglement between 
opposite momenta, as expected from a superconductor. 

It is however, more interesting to study the entropy 
associated with the reduced density matrix p' as one ap¬ 
proaches a topological transition. In Fig. [^we consider 
the dependence of the entropy with momentum for differ¬ 
ent values of the chemical potential, ep = —3,= 1, for 


different values of the magnetization, Mz^ as one crosses 
the transition lines between a C = 0 phase and a C = 1 
phase and between a C = —2 and a C = — 1 phases. As 
with other signatures the time-reversal momentum val¬ 
ues play a special role associated with the transition lines. 
In the top panel of Fig. pT] there is a change of behavior 
around momentum (0,0) and in the lower panel around 
momentum (7r,7r). From the data analysis, we can say 
that in the (7 = 0 phase the entropy has a small, fixed 
but finite value, at zero momentum. As one crosses the 
transition to the (7 = 1 phase the entropy vanishes at 
this momentum. Also, as one approaches the transition 
line the entropy slope increases around zero momentum. 
After the transition the slope decreases considerably and 
the entropy is smaller than in the phase with zero Chern 
number. The lower panel shows a similar behavior but 
around momentum (7r,7r), the slope becomes quite large 
at the transition and then decreases. Before the transi¬ 
tion, as the magnetization is increasing, the entropy is 
also finite (but rather small) and after the transition it 
vanishes at that momentum. Also, the magnitude of the 
entropy is smaller after the transition, since the magnetic 
field is larger. 

In Appendix B we present results for the eigenvector of 
the largest eigenvalue of the /c-subspace reduced density 
matrix. 



4 0 


FIG. 12: (Color online) Total fidelity as a function of chem¬ 
ical potential and magnetization. The fidelity was calculated 
using two points from parameter space, Mz and Mz H- SMz, 
with 5Mz = 0.01. Note that the vertical scale is inverted. 


V. FIDELITY AND FIDELITY SPECTRUM 

The quantum fidelity between two pure states (for two 
sets of parameters) is the absolute value of the overlap 
between the ground states for the two sets of param¬ 
eters. In general, the quantum fidelit}® between two 
states characterized by two density matrices pi and p 2 ^ 
may be defined as the trace of the fidelity operator, 7^, 
as F{pi,P 2 )= TrJ" = 

The fidelity operator T can be studied using differ¬ 
ent basis states, associated with different representations, 
such as position, momentum, energy or charge and spin. 
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Rewriting the fidelity operator in these different repre¬ 
sentations, allows us to look more directly at the spe¬ 
cific relevant modes that participate more actively in the 
critical phenomena accompanying the phase transition. 
In this way, we can obtain a more complete and phys¬ 
ical characterization of the phase transition and of its 
underlying physics mechanisms. While the entanglement 
spectrum has some relation to the energy spectrum of the 
edge states or even bulk states, the fidelity spectrum con¬ 
tains information about which eigenvalues have a larger 
contribution to the distinguishability between quantum 
state^^. 

Considering two density matrices, that result from the 
momentum space partition, for two points in parameter 
space, we can write that 

pikQik = ; p2kQ2k = Q2k^2k (23) 

As mentioned above, the fidelity between two states, 
characterized by two density matrices and p 2 k^ iiiay 
then be defined as the trace of the fidelity operator, Jbc, 

Fk{pik^P2k) = FvTk = ( 24 ) 

The square root of the density matrix can be written as 

^k = Qik^kQiP (25) 

and 

j^k^ = \/p^p2kV^ ( 26 ) 

Diagonalizing 

Tk = UkfkUp (27) 

we write that 

j^k,n = Uk,nfk,nUj^l^ ; n = 1,..., 16 (28) 

Therefore the /c-fidelity is obtained as 

Fk=TT{J^k) = J2fk,n (29) 

n 

and the total fidelity is finally obtained as 

F = llFk ( 30 ) 

k 

We begin by considering that the two density matrices, 
Pik^ P2k^ points nearby in parameter space. 

As the fidelity measures the distinguishability between 
different states, as one approaches a transition one of the 
density matrices is calculated on one side of transition 
and the other is either at the transition point or in the 
other phase (depending on the step in parameter space) 
and the fidelity shows a minimum. 

In Fig. [^we present results for the total fidelity as a 
function of the chemical potential and the magnetization. 
Far from the transition points the fidelity is close to one. 



FIG. 13: (Color online) Highest eigenvalue of the /c-fidelity 
spectrum at a transition point sf = —1, Mz = 1 as a function 
of momentum. The fidelity was calculated using two points 
from parameter space, sf and sf + ^Sf^ with Ssf = 0.01. 
The transition occurs for momentum k = (tt, 0) and equiv¬ 
alent points, as shown by the decrease of the fidelity in the 
neighborhood of the singular momenta values. 


as the two states described by the two density matrices 
are very similar. The minima in the fidelity faithfully 
track the transitions previously described (note that the 
vertical axis is upside down). As for the entropy, fixing 
the magnetization at = 0.1 the /c-fidelity faithfully 
singles out the singular momenta values. 

The transitions may also be signaled by the decrease 
of the highest eigenvalue of the /c-fidelity spectrum at 
the transition point, when the two density matrices are 
calculated at points in parameter space surrounding the 
transition. This is shown in Fig. for a particular 
example, Sf = —1, Mz = 1, for which the transition line 
is associated with the singular points k = (7r,0), and 
equivalent points in the Brillouin zone. 

Complementary, interesting information, can be ob¬ 
tained taking the two density matrices in the fidelity ex¬ 
pression, corresponding to points in parameter space that 
are far apart. This allows to compare different phases, 
and not specifically to detect the locations of the phase 
transitions. It provides interesting information about the 
nature and momenta values responsible for the transi¬ 
tions. Also, it describes the distinguishability between 
two phases in momentum space. 

In Figs. |14|15| we present results for the /^-fidelity op¬ 
erator spectrum, where the two density matrices define 
two states in parameter space that are deep inside dif¬ 
ferent topological phases. In Fig. we consider one 
density matrix at a state specified by a point in a trivial 
phase with (7 = 0 (specifically ep = —l^Mz = 0.5) and 
the other density matrix corresponds to states at various 
magnetizations, along the vertical dashed line of Fig. 
for which the Chern numbers are C = —2,(7 = —1,(7 = 0 
and the magnetization is Mz = 2, Mz = 4, Mz = 5.5, 
respectively. As we move along that line in parameter 
space, we cross, in sequential order, lines characterized by 
momenta k = (7r,0),A; = (0,0),/c = (7r,7r), respectively. 
The /c-space fidelity highlights the differences between 
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FIG. 14: (Color online) /c-fidelity operator spectrum Fk for 
PA corresponding to = —1, Mz = 0.5 where the Chern C = 
0, and pb corresponding to £f = — 1 as well and Mz — 2^ C — 
—2, Mz = 4, C = —1 and Mz — 5.5, C = 0, respectively. 


the various phases. Its deviations from unity throughout 
the Brillouin zone are significant, even when the phases 
have the same Chern number, as a consequence of the 
different band-fillings. Singular transition points trans¬ 
late to zeros in the /^-fidelity operator spectrum. As the 
top panel shows, the /c-fidelity operator spectrum van¬ 
ishes at the momenta k = (tt, 0) and equivalent points 
(two of them are independent). The middle panel shows 
in addition a zero at the center of the Brillouin zone 
k = (0,0). Consequently the last panel shows zeros at 
all singular momentum locations. The /c-fidelity opera¬ 
tor therefore signals the minimum number of transitions 
that must occur when going from one trivial phase to the 
other trivial phase. 

Similar results are shown in Fig. In this case 

we also start from the same trivial phase with (7 = 0 
{sf = —2^Mz = 0.5) and consider a sequence of den¬ 
sity matrices, associated with states along the horizon¬ 
tal dashed line of Fig. In this instance we consider 
in the top panel that the second density matrix corre- 


FIG. 15: (Color online) /c-hdelity operator spectrum Fk for 
pA corresponding to £:f = —2,Mz = 0.5 where the Chern 
(7 = 0, and ps corresponding to Mz — 0.5 as well and 
sf = 0,(7 = —2, £f = 2, (7 = 0 and at bottom panel pA 
corresponding to £f = 2,Mz =0.5 where the Chern (7 = 0 
and pB corresponding to Mz = 0.5, stf = 6,(7 = 0, respec¬ 
tively. 


spends to Mz = 0.5, = 0, (7 = —2. A line asso¬ 

ciated with k = (7r,0) is crossed. In the middle panel 
we consider that the second density matrix is at point 
Mz = 0.5, = 2, C = 0. A second transition line at the 

same momentum is now crossed, which implies a zero of 
quadratic order, since the same momentum becomes gap¬ 
less twice. A similar result occurs if we consider a density 
matrix at Mz = 0.5, = 2 and the other at the same 

magnetization and = 6. In this case the momentum 
k = ( tt , tt ) contributes twice with a gapless spectrum and 
the dispersion is now quadratic around that momentum 
(the corners of the Brillouin zone), as shown in the bot¬ 
tom panel of this figure. 
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VI. CONCLUSIONS 

Various quantum information techniques have been 
used in this work to detect and distinguish topological 
phase transitions in superconductors. In particular, the 
entanglement von Neumann entropy, entanglement spec¬ 
trum, fidelity and the recently introduced fidelity spec¬ 
trum. In addition to a faithful identification of the tran¬ 
sitions, these methods clearly signal the relevant modes 
responsible for the transitions. The analysis of the high¬ 
est weight eigenvector components clarifies the nature 
of each phase and the contents associated with a tran¬ 
sition between two phases. Also, allows to distinguish 
between a trivial and a topological phase and between 
the Z 2 and Z topological phases. Specifically, for the 
p-wave superconductor considered here, the spin triplet 
basis state has a particular role in the non-trivial topo¬ 
logical phases reaching saturation. A detailed analysis of 
the momentum dependence of the coefficients of the ba¬ 
sis states of the highest state eigenvector was carried out. 
The analysis of the /c-fidelity operator spectrum allows a 
characterization of the sequence of the minimal number 
of transitions that have to occur between two arbitrary 
points in the phase diagram. 
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Appendix A: Highest weight density matrix 
eigenvector 

In this appendix we study in detail the coefficients of 
the highest weight eigenvector. We consider a sequence of 
states as a function of the chemical potential, for dif¬ 
ferent values of the magnetization Mz. We compare the 
results for Mz = 0 and Mz = I for ep = —4,—2,0,2,4, 
typically. Along these sets of values the Chern number 
changes from topological phases with (7 = 0 to phases 
with finite Chern number. As discussed in the main 
text, even though some phases have no Chern number, 
there are edge states except in the (7 = 0 phases with 
small magnetization and \sp\ > 4, and large magnetiza¬ 
tion around sp = 0^ indicated in Fig. For instance, 
considering Mz = I and changing the chemical poten¬ 
tial, the Chern goes through a sequence of values, C = 
1,(7 = 0,(7 = —2, (7 = 0, (7 = I, from non-trivial to triv¬ 
ial phases, crossing in sequence lines that are singular at 


the momentum values k = (0, 0), /c = (tt, 0), /c = (tt, tt). 

In Fig. [^we present, in a rather compact way, the ab¬ 
solute values of the coefficients of the 16 basis states, for 
the highest weight eigenvector, as a function of momen¬ 
tum, ordered sequentially in the two-dimensional plane. 
We compare the results for Mz = 0 and Mz = I for 
different chemical potentials. The color codes of the co¬ 
efficients is the same as in Fig. Clear differences can 
be seen between the Z 2 phases and the Z phases. For 
Sp = —4 and Mz = I the system is in a phase with 
(7 = 1 but for Mz = 0 the system is at the frontier be¬ 
tween a completely trivial phase with (7 = 0 and no edge 
states and a Z 2 phase also with (7 = 0 but edge states. 
The dominant state is the state |0000) in most of momen¬ 
tum space. For Mz = 0 the two spin orientations of the 
triplet state, |II00) and |00II), are degenerate and are 
superimposed. The other state that gives a significant 
contribution is the fully occupied state |IIII). 

For values of the chemical potential, \ep\ > 4, there 
is a clear gap between the dominant empty state and 
all others that is particularly large when Mz = I (these 
results are not shown). At the transition to the Z 2 phase 
the gap closes at some momentum values. The absolute 
value of both triplet components saturates at 1/2. This 
is to be contrasted with the case when there is a finite 
magnetization. First, in this case the degeneracy between 
the two triplet states is lifted. Inside the phase with finite 
Chern number, the triplet component |II00) saturates to 
one at the appropriate singular momentum value (zero 
momentum in the case of ep = —4). This saturation 
is characteristic of the phases with finite Chern number. 
Increasing the chemical potential, for instance Sp = —2, 
the weight of the state |IIII) increases and the weight 
of the spin triplet state decreases. For Mz = 0 there is a 
slight decrease but a clear difference is noted when Mz = 
I. In the topological phase the triplet state is dominant 
at some momentum values while in the trivial phases 
(with (7 = 0) the contribution from other states is also 
important. At ep = 0 {C = —2) there is also a saturation 
of the triplet component, but at different momenta. It 
is also clear that, as the chemical potential increases, 
the weight of the fully occupied state increases, until it 
becomes dominant over the Brillouin zone. The roles of 
the states |0000) and |IIII), are naturally reversed as the 
chemical potential changes form the bottom to the top 
of the tight-binding band. 

To further understand the relative weights of the basis 
states we show in Figs. the momentum depen¬ 

dence of the sum of the absolute squares of i) empty and 
doubly occupied states, ii) singlet pairing states and iii) 
triplet pairing states. As mentioned before, these are all 
the states that give significant contributions. This allows 
to better compare the cases Mz = 0 and Mz = I be¬ 
cause the spin orientations become degenerate. (Also we 
stress that previous results were presented for the abso¬ 
lute value of the coefficients, while here we consider the 
absolute values squared). In Fig. 17 we consider Mz = 0 
and in Fig, we consider Mz = I. The various rows 
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FIG. 16: (Color online) Absolute value of the coefficients of the basis states of the highest weight eigenvector as a function of 
momenta values. In the left panels Mz — 0 and in the right panels Mz — 1. The rows are for the values of = —4, —2, 0, 2, 
respectively, from top to bottom. The color codes of the coefficients is the same as in Fig. 


correspond to different chemical potentials. In the left 
column we plot the sum of the (square of the) empty and 
doubly occupied states, in the middle column the singlet 
pairing and in the right column the triplet pairing states. 

The results shown in these figures summarize the var¬ 
ious aspects found previously and further clarify new as¬ 
pects. a) The empty and doubly occupied states usually 
dominate over the Brillouin zone, for all parameter val¬ 
ues, except near the Fermi surface, b) Near the Fermi 


surface the main contribution comes from the triplet and 
singlet pairing states, c) Since we choose the amplitude 
of the s-wave gap parameter smaller than the triplet am¬ 
plitude, the singlet pairing coefficient is smaller than the 
triplet. For Mz = 0 the amplitudes of the singlet and 
triplet states are similar. For Mz = 1 the singlet states 
contribution is considerably smaller than the triplet one. 
d) In the case Mz = 0 the maximal amplitude of both the 
singlet and triplet states is quite uniform along the Fermi 
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FIG. 17: (Color online) Sum of the absolute square values of the coefficients for the empty state and doubly occupied state 
(left panels), for the spin singlet states (middle panels) and triplet states (right panels), for Mz = 0 and for sf = —4, —2, 0, 2,4 
from top to bottom, as a function of momentum. 


surface with small peaks near the singular momenta, e) 
In the case of finite magnetization, the triplet coefficient 
saturates near the singular points. Therefore, the depen¬ 
dence on momentum along the Fermi surface is stronger, 
f) As for the entanglement measures, the significance of 
the singular momenta is quite clear in the various topo¬ 
logical phases. 


Appendix B: Eigenvectors of /c-subspace reduced 
density matrix 

The /c-subspace reduced density matrix p' is defined in 
a basis of the occupation numbers The states 

are therefore of the type |00), |01), |10), |11). The diag- 
onalization of this 4x4 matrix leads to the eigenvalues 
and eigenvectors expressed in this basis. 

In the first two panels of Fig. we present results 
for the eigenvector of the largest eigenvalue for the point 
in parameter space Mz = 2,5^ = 0 (C = —2). We plot 
the sum of the absolute coefficients along the states |00) 
and 111), in the top panel and along the states |01) and 
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FIG. 18: (Color online) Sum of the absolute square values of the coefficients for the empty state and doubly occupied state 
(left panels), for the spin singlet states (middle panels) and triplet states (right panels), for Mz = 1 and for sf = —4, —2, 0, 2,4 
from top to bottom, as a function of momentum. 


110), in the lower panel. We see that the contributions 
from the states with empty or double occupancy and the 
states with single occupancy, are exclusive which implies 
a very sharp separation in the Brillouin zone between the 
two sets of states. This rather sharp separation occurs 
around the phase diagram. 

However, a detailed analysis near a transition between 
the various topological phases shows that, focusing on 
the special moments involved in a given transition these 
provide a clear signature of the crossing. As an exam¬ 
ple we consider in the last three panels of Fig. the 
eigenvector components along the states |00) and |11) of 


the highest eigenvalue of the reduced density matrix for 
= — 3 and Mz = 0.9,1.0,1.1, respectively. This tran¬ 
sition is associated with the gapless point at momentum 
(0,0). As Fig. [l^ shows there is a clear signature at 
this momentum value as the transition occurs. In the 
regime when (7 = 0 there is a peak of the eigenvector 
components at k = (0,0) which vanishes in the phase 
with (7 = 1, a typical behavior of the time-reversal mo¬ 
menta at the transition points. We have checked that in 
the components of the singly occupied states the same 
(complementar) feature occurs. 
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FIG. 19: (Color online) First two panels: Eigenvector components of the highest eigenvalue of the reduced density matrix for 
Mz = = 0 (C = —2). Following panels: Eigenvector components along the states |00) and |11) of the highest eigenvalue 

of the reduced density matrix for sp = —3 and Mz = 0.9,1.0,1.1, respectively. 
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